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ABSTRACT 


The development of new decision support systems for 
Antisubmarine warfare will entail the installation of 
propagation loss models on ASW aircraft. The decision to put- 
either a range dependent or range independent model in the 
system will affect the predicted ranges, the overall 
probability of detection, and the computation time. 

Comparisons of the range dependent ASTRAL and range 
independent RAYMODE propagation loss models were made in the 
Eastern Mediterranean, the Gulf of Oman and the South China 
Sea for eight source/receiver/frequency combinations. 
Computation time differences between the two models were not 
Significant at either of the source frequencies (50 Hz or 400 
Hz). RAYMODE showed much better correlation with the split 
step PE model which was used as a standard. The ASTRAL model 
often predicted lower transmission losses than either RAYMODE 
or PE. For the short detection ranges normally encountered in 
air ASW the more complex range dependent models are not 
necessary. The RAYMODE model or a comparable range 


independent model will provide adequate propagation loss 


predictions. 
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I. INTRODUCTION 

The intent of this paper is to present a comparison study 
of the effect of using a range dependent versus a range 
independent propagation ioss model in an airborne ASW 
tactical decision support system. In an ASW tactical 
Situation, processing time and accurate probability of 
detection are extremely important. The ability to produce 
quick and reliable transmission loss information while on 
station would be a definite asset in a decision support role. 
The two transmission loss models currently being considered as 
candidates to be installed in an aircraft decision support 
system are the range dependent ASTRAL or the range independent 
RAYMODE model [Ref. 1]. Using a range dependent (RD) versus 
a range independent (RI) propagation loss model may affect 
both the processing time and the probability of detection. 

This study compares the ASTRAL and RAYMODE models in 
several areas. First, it will review the physics that both 
models use to develop their propagation loss information. 
Next, using three different geographical areas and a 
combination of two source frequencies, two source depths, and 
two hydrophone depths, transmission loss profiles will be 
generated. The profiles are then used to develop detection 


ranges and probabilities of detection for three specific 


sonobuoy patterns (2x8, 4x4, 5x6x5, shown in Fig. 1, spacing 





is based on the predicted detection ranges). The comparisons 
of the produced transmission loss profiles, the detection 
ranges, and the probabilities of detection of the patterns 
will all be presented. The processing time required by each 


model to run the transmission loss profiles will also be 


discussed. 





- Figure 1 - Sonobuoy Pattern Geometries 





II. PROPAGATION LOSS MODELS 

The Navy currently uses many different oceanographic and 
propagation loss models. Propagation loss is the attenuation 
of the sound intensity as it passes through the ocean from the 
source to the receiver. Propagation losses are, in general, 
affected by such factors as the temperature/depth profile of 
the water column from the surface to the bottom, the source 
frequency, the condition of the ocean surface (wind speed or 
sea state), bottom depth, bottom composition and thickness, 
salinity, source and receiver depths, and range to the source. 
Once a model is tested and accepted by the Navy it is 
designated as a Navy Standard Model. There are both range 
dependent and range independent standard models . The range 
independent models assume the water column does not change 
with range and do not take into account the ocean bottom 
topography. A single sound speed profile, water depth and 
bottom loss are used for the calculations for an entire area. 
The range dependent models do use bottom eopogeaphy and allow 
the sound speed profile and bottom loss to change with range. 
RAYMODE is currently the only range independent Navy Standard 
model, while ASTRAL and PE (Parabolic Equation model) are the 
most commonly used RD models. The RAYMODE and ASTRAL models 


use distinctly different methods for calculating the 


transmission loss from a source to a receiver. In general, 











most propagation loss models rely on either wave theory or ray 
theory (Ref. 2:p. 1]. 

. The wave equation is a partial differential equation 
relating the acoustic pressure p to the coordinates x,y,z and 


the time t, and may be written as 


B® . ga, BD, Bp) : 
ez So Bye *By2 Gaz? = 


where c is the sound speed and may vary with the coordinates. 
There are two basic theoretical approaches to a solution of 
the wave equation for the specified boundary conditions. The 
first is wave theory, in which.the propagation is described in 
terms of characteristic functions called normal modes, each of 
which is a solution of the wave equation. The normal modes 
are combined additively to satisfy the source conditions of 
interest. The result is a complicated mathematical function 
which, though adequate for computations on a computer, by 
itself gives little insight on the distribution of the energy 
of the source in space. The other basic approach to a 
solution of the wave equation is ray theory. The essence of 
ray theory is that the propagation path of individual rays can 
be traced by using Snells Law. The reader interested in an 
introductory description of wave and ray theory is referred to 
a book by Officer [Ref. 3]. Snell’s Law, one of the most 
important practical tools of ray theory, describes the 
refraction of sound rays in a medium of variable sound speed. 


As applied to an ocean of horizontal layers of isospeed water, 





Snell’s law states that the grazing angles 0,, 0,,...of a ray 
at the layer boundaries are related to the sound speeds c,, 


Co,..- Of the layers by 


cos cos® cos 
t= 22 = — 3 = ,...=a constant for any one ray 


Cy C2 C3 





This expression enables a particular ray to be traced out by 
following it through the successive layers into which the 


sound speed profile may have been divided. [Ref. 4] 


A. RAYMODE 

The passive RAYMODE model predicts the transmission loss 
of sound through an ocean medium. It is a Navy Standard 
Passive Acoustic, Range Independent, Propagation Loss Model. 
The model was developed at the Naval Underwater System Center, 
New London, by Dr. G.A. Leibiger in 1979 [Ref. 5:p. 1]. 
RAYMODE has undergone many revisions since it was first 
presented, but little has been written about the model since 
1985. The Navy is currently using Revision 8.1 of the RAYMODE 
program [Ref. 6:p. iv]. Although the program has been updated 
to simplify user interface and to utilize new computer 
technology, the basic physics of the model for calculating 
propagation losses has changed very little. The RAYMODE model 
combines the wave and ray theories, described above, using a 


normal mode summation technique for each propagation path. 


RAYMODE gives control of many of the propagation loss 





variables to the user to provide the most precise definition 
of the environment. The propagation losses are computed for 
each path type (e.g., direct path, convergence zone, or bottom 
bounce), at user-defined ranges and summed for each source 
angle, or for each of a group of normal modes. The total 
field at each range is determined by summing the contributions 
over all path ivose. 

| The operation of the RAYMODE model as described by R. C. 
Medeiros in his overview of the model [Ref. 2:pp. 4-12] is 
summarized below: 

1. The operator is required to provide environmental data 
and program control parameters, or indicate which historical 
Gatabase to use. These inputs include: 

a.__Environmental data: 
Sound speed versus ocean depth profile 
Bottom depth 


Bottom type code (BLUG) 
Wind speed 


b. _Geometrical data: 
Source depth 
Receiver depth 


c. Source frequency 


dad. Performance data: 
Minimum and maximum performance ranges 
Performance range step 


2. Input data are processed by the program in preparation 
for the actual computation and summation of propagation 


losses. 





a. Propagation sound speeds are computed which 
correspond to the input source and receiver depths by linear 
interpolation between the entered profile depths. These 
depths and their associated sound speeds are inserted into the 
warding profile. Together, they define the end-points of the 
propagation paths. 

b. The bottom depth is inserted into the waevSng 
profile and a corresponding sound speed is extrapolated or 
interpolated. The profile is then modified by the program to 
account for the effect of earth curvature. 

c. A gradient profile is computed from adjacent values 
of c as a function of depth. The sound speed gradient, 
g = Ac/Az, is the ratio of the change in sound speed and the 
depth change over which the sound speed change occurs. 

ad. The profile is thus partitioned into wave number 
domains, (see Fig. 2) [Ref. 2:p. 5] testing the profile for 
localized maxima and minima which indicate those water column 
characteristics, i.e., ducts and layers, that define the 


propagation paths. The wavenumber 


wavenumber = —Frequency(rad/sec) | (3) 


soundspeed (yard/sec) 


is the parameter which most effectively specifies propagation 


path characteristics such as propagation angles and depths and 


ranges associated with different propagation modes. 





ome 


—*— DEPTH (yds) 


Sound Speed (Y¥d/sec) ——-———>= 


Sound speed at the surface 
Sound speed at the source depth 
Source depth 
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Local maximum 

Local maximum 

Bottom depth 
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Figure 2 - Sound Speed Profile Indicating Criteria for 
Wavenumber Partitioning 








3. Once these wavenumber domains are determined, they are 
interpreted by the program as angular intervals, a bounded set 
of sound propagation angles. 

4. For each path, two more quantities are defined. The 
first is the number of cycles that exist along the path range 
between the source and receiver. A full cycle is defined as 
the path from surface or apex to the bottom or nadir and 
return. A partial cycle is the path from receiver to surface 
or from the source to surface. The total number of cycles and 
the partial cycles when appropriately combined yield the path 
range for a given propagation path. Also a phase change 
occurs at the turning points which must be accounted for. One 
additional parameter to be computed is the number of modes 
which exist for the wavenumber interval under examination. 
For some paths there may thousands of modes, assuming the 
frequency is sufficiently high. 

When the number of modes propagated in an interval 
becomes too large, the interval is subdivided into what are 
called K-regions. A K-region is the interval between two 
wavenumbers that defines a propagation path. Referring back 
to Fig. 2, the sound speed profile is partitioned into 
intervals bounded by a minimum and maximum sound speed. By 
converting these sound speeds into wavenumbers the interval 


between the smallest and largest wavenumber is subdivided into 


a predetermined number of K-region intervals. The K-region 





interval is then converted into a propagation angle using 
Snell’s law. 

Each K-region interval corresponds to an initial angle of 
a propagated ray that is extended down the range path with 
curvature based on the application of Snell’s law to the sound 
speed profile. At each range increment the sum of all the 
mode contributions from each K-region interval is made 
resulting in a single combined value of acoustic pressure 
field. The individual contributions are summed at each range 
for all the rays in an angular interval, which is then summed 
to those of previous angular intervals to obtain the total 
acoustic pressure field at a given range. After correction 
for boundary losses and sonar beam deviation loss the result 
is converted to decibels and output as propagation loss 
versus range. [Ref. 1l:pp. 4-6] 

When the propagation angle increases, as it does in 
bottom bounce regions, the number of modes propagated rises, 
increasing the number of ieeeations necessary. To avoid this, 
when the number of modes exceeds 10, the model integrates over 
a series of K-regions and power sums or coherently sums the 
contributions of each. 

Two types of propagation losses are computed, coherent 
and random incoherent. Coherent losses are those attributed 
to ray paths with an exact phase relationship between all 
paths. Random losses assume a random distribution of phase. 


Because coherent propagation losses tend to vary widely, small 


10 








modifications in input data can cause large changes in output 
data. The operational RAYMODE model utilizes the random phase 
propagation loss predictions outputs which tend to be 
smoothed over the performance range span, providing an 


expected value of propagation loss. [Ref. 2:p. 12] 


B. ASTRAL 

ASTRAL (ASEPS Transmission Loss) is a range dependent, 
low frequency propagation loss model developed by Seienee 
Applications International Corporation [Ref. 7:p.ii}. Since 
ASTRAL’S original version in 1979 it has undergone several 
modifications. The addition of a surface duct and convergence 
zone capability has increased its applicability and utility. 
ASTRAL 4.1 is currently used as one of the Navy Standard Range 
Dependent Propagation Loss models [Ref. 7:p.3}. The remainder 
of this section is a summary of the basic physics and theory 
of the ASTRAL model as described in Refs. 7 and 8. 

The calculations in ASTRAL are carried out in two parts, 
the near field and the far field. The near field calculations 
are based on a ray trace scheme. The number of rays traced is 
set to 50 in the operational version. The rays are traced out 
to a range specified by the user. The angular spread between 
rays is adjusted so that 80% of the rays have launch angles 


less than 45 degrees and 20% have launch angles greater than 


45 degrees. 





In the far field each ray is associated with a set of 
modes or "“smodes". These smodes approximate the envelope of 
the oscillatory depth function of true normal modes, where 


pressure can be expressed in the form 


P = A,%,U,(Z,) U,(Z,) exp (1k,r) / /K,F (4) 
where %, and 4% are the source and receiver depths 
respectively. The basic ASTRAL model begins with the normal 


mode equation [Ref. 7:p. 12] 


P(e) 2) SS Kau (2) OZ.) exp (ik,r) (5) 
where S, is a constant related to the source strength divided 
by the square root of the range r, the U are normalized depth 
functions, k, is a reference wave number (equal to 2nt/co , 
where c, is a reference sound speed), and the asymptotic form 
has been used for the Hankel function. The symbol i is the 
square root of minus one. The intensity is assumed to be in 
the form 

PSA VAs Vee 7K; (6) 
where 
y2(z) = [z,(6,) tan6,(z)J~ (7) 
and where r,(0,) is the apex-to-apex cycle range of a ray 


launched with angle 6,. The summation over modes is replaced 


by a summation over sets of modes, so Equation 4 can be 


written as 





TAD ee.) Yee) x, (8) 
where each value of the running index corresponds to a set of 
modes, or "smodes." The variable X, is given by N,/k,, where 
N, is the number of modes in the smode. In the ASTRAL code, 
these average values are actually calculated at specific ray 
angles measured relative to the axis of minimum sound speed 
and thus a smode corresponds to an angle bin. Using the piiaide 


integral ¥,, where 


Vn, = ftsin(®,(2)) /(2z) } dz (9) 


then N,=0,,/". If w,/(sin63@) is assumed constant, then the 


pressure can be expressed as 


DB = Dam 25) Vn Ze) Xqgko eM? (10) 


where 
Yn(Z) ~ f¥z(z) (11) 
and 


X, = cos[o,¥,,(2,) ]cos[w.p,,(Z,)]exp (ik,r) . (12) 


The smode is then extended at the turning points via Airy 


functions. The smodes are propagated in the far field until 


an environmental change occurs. If the structure of the 





sound-speed field changes then a new smode assignment must be 
made using the new turning point structure and assuming 
adiabatic invariance. As the environment changes, the phase 
integral is kept constant by choosing a new turning point 
sound speed. [(Ref. 8] 

Convergence zone effects can be added if each smode is 
assumed to be composed of four parts that correspond to the 
four different up/down arrivals at both source and receiver. 
Each of the components is propagated in range using the cycle 
length of the smode and the intensity of each is assumed to be 
proportional to the square of a sin(x)/x function. At each 
range the intensity is summed and smoothed.[Ref. 8] The full 
derivation of this concept is beyond the scope of this paper. 


For a full description of Convergence Zone ASTRAL the reader 


is referred to Ref. 7, pp. 14-18. 








III. BACKGROUND 


A. DESCRIPTION OF AREAS USED 

Three areas were chosen for the ASTRAL and RAYMODE 
comparisons. The areas were chosen based on the new 
regionally focused Navy strategy stated in "From the Sea...". 
The specific coordinates for the areas investigated were 
picked at random and then rounded to the nearest degree for 
ease of data entry. Each of the three areas are 50 nm by 50 
nm boxes. 

The first area is in the Eastern Mediterranean centered 
at 33N/033E. Area 1, located between Egypt and Cyprus as 
shown in Fig. 3, is typical of one in which an air ASW 
platform might have to search for a third world eicaarine: 
The area is close enough to the coast to be sonsiaeed 
littoral but too deep to really be considered shallow water. 
The bottom contour and sound speed profile are shown in Fig 4. 
The surface temperature in the area is 60° F , there is a 
layer depth of 350 ft, the depth excess is 525 fms and there 
is a sound channel at 715 ft. 

Area 2 in the Gulf of Oman, is centered on 24-30N/059E. 
This region has become a fairly high interest area since the 


Persian Gulf War. The specific area used for the study is 


shown in Fig. 5. The sound speed profile for the center of 








Figure 3 - Location of Area 1 in Eastern 
Mediterranean 


Bottom Profile of Area i 
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Figure 4 - Sound Speed Profile and Bottom Profile for Area 1 
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Sound Speed Profile for Area 2 
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Figure 6 - Sound Speed Profile and Bottom Profile for Area 2 
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Figure 8 - Sound Speed Profile and Bottom Profile for Area 3 
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the area, and the bottom profile are shown in Fig. 6. A 
surface temperature of 80° F results in a large negative 
gradient below the layer depth of 250 ft. There is no depth 
excess (depth needed for convergence zone propagation to 
exist) in this area and the sound channel is fairly deep at 
approximately 4000 ft. . 

The third area is in the South China Sea at 10N/110E, 
only 75 miles off the coast of Vietnam (Fig. 7). The sound 
speed profile is similar to the one in Area 2 with a strong 
negative gradient below the layer (Fig. 8). The surface 
temperature is 75° F, and the sound channel is at 2000 ft. 


Similar to Area 2, there is no depth excess. 


B. DATA COLLECTION 

The ASTRAL, RAYMODE, and PE models used for this study 
were run inside a new developmental ASW Tactical Decision Aid 
(ASWTDA version 2.1.2.1) which has many ASW tactical 
applications. The ASWTDA program is run on a SUN work station. 
The databases used by ASWTDA are the Oceanographic and 
Atmospheric Master Library (OAML) and the GEM database. The 
GEM database is a subset of OAML and was created by ASWTDA‘’s 
developing contractor, SONALYSTS, Inc., in order to reduce the 
OAML database size and the time required to access the 
information. Table I lists the versions of the models and 


databases used by ASWTDA. [Ref. 9] 
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ASWTDA (V.2.1.2.1) } 

















| 
ESS <(-" ee e ( ee 
|___SHIPNOISE | 4 


PASSIVE RAYMODE 
Pee oe SO 
CZ ASTRAL 


Table | - Models and Databases Used by ASWTDA 





LFBL and HFBL are low and high frequency bottom loss 
databases, GDEM is an oceanographic environmental database and 
DBDBC is the bathymetry database. The areas of interest were 
generated inside ASWTDA's “Environmental” function with 20 
sample points inside each 50 nm by 50 nm box. Propagation 
loss profiles were generated for the combinations of 90 and 
400 ft receivers, 150 and 450 ft depth targets, and 50 and 400 
Hz sources. This produced 8 range dependent and 8 range 
independent profiles for each area. The specific technique 
and code used to extract this data is shown in Section B of 
Appendix A. The propagation loss profiles were then evaluated 
using a Matlab subroutine to produce Median Detection 
Ranges(MDR) for Figures of Merit (FOM's) of 65,70,75,and 80 
adB. Three distributive sonobuoy patterns (4x4, 5x6x5, 2x8) 


were evaluated at each FOM and receiver/source/ frequency 
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combination. A pattern evaluation program provided by Wagner 
Associates called NESST was used to produce probabilities of 
detection for the three selected patterns. A description of 
the NESST program and how the CDP's were obtained is included 


in Appendix B. 


Cc. ASSUMPTIONS 

For the purpose of this study several basic assumptions 
were made. The first is that the split-step PE propagation 
loss model is sufficiently accurate to serve as the basis for 
comparisons. The target, for the probability of detection 
calculations, is assumed to be uniformly distributed in a 50 
nm by 50 nm box. In most littoral ASW situations the passive 
detection ranges are relatively short (1 to 10 nm). Therefore, 
even though the data for the entire area was used for analysis 
purposes, the profiles generated were only graphed out to 30 
nm. The sonobuoys were considered to be cookie cutter 
detectors: if the target was inside the given range of the 


specific sonobuoy it was counted as a detection and if it was 


outside the range it was considered as undetected. 





IV. ANALYSIS 


A. AREA i: EASTERN MEDITERRANEAN | 

Figures 9-12 compare the ASTRAL and RAYMODE propagation 
loss profiles against PE generated profiles. Analysis of the 
ASTRAL profiles shows much lower dB loss versus range than in’ 
the PE model or the losses expected from spherical spreading. 
(The "x" marks on the Siache indicate spherical spreading at 
1,2, and 3 nm.) The RAYMODE profiles correlate very well with 
the spherical spreading and the PE profiles, with the 
exception of the 90ft/150ft/400 Hz run. This lack of 
correlation in the shallow/shallow situation is most likely 
due to the fact that both the source and receiver are in the 
mixed layer and the 400 Hz source frequency is above the 150 
Hz cutoff frequency for a layer of that depth. The RAYMODE 
model aces not predict the surface duct and chevetore gives 
much shorter predicted ranges. 

Predicted detection ranges were obtained from the 
proploss neseiiee by extracting the range at which the profile 
crossed the given figure of merit. The ranges for each FOM 
are tabulated in Appendix A. The average differences in the 
detection ranges are shown in Table II. 

The calculated detection ranges were used to construct 


the three sonobuoy patterns, with a buoy spacing of 1.5 MDR 
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___ ASTRAL - PE 
| Run} FOM6s | 
LAIMMH [02 
PAIMML [0.3 
TAIMSH [70.2 
PAIMSL [0.4 
PAISMH | 0.4” 
PAIsML [0.5 
PAISSH [0.2 
PAisst [0.6 
AVG DIF ; 


| Run | FOM70__| 
[AIMMH |" 0.8 
AIMML [1.8 
PAIMSH [0.5 
PAiMSL [14 
PAISMH [1.0 
PAiSML [1.4 
PAiSSH | 


| 02 
A1SSL 
AVG DIF 1.0 


| Run_ | FOM75_| 
PAIMMH | 3.0 
PAIMML | 3.9 
|AIMSH | 2.0 
PAIMSL | 2.9 
}AISMH | 2.6 
PAISML | 4.9 
PAISSH | -16 
AVG DIF 


Differences 


RAYMODE - PE 


FOM 70 


pS 
a <a 
07 
4.2 


FOM 75 
a ee 
0.2 
0.7 


ASTRAL - RAYMODE 





Table II - Area 1 Average Predicted Detection Range 








and a row spacing of 1.67 MDR as prescribed in the ASW 
tactical manuals. The pattern in each scenario that generated 
the largest cumulative detection probability (CDP) was used in 
Fig. 13 to compare the CDP's of each model. (The first column 
in the table is the run ID which is described in more detail 
in Appendix A.) It should be noted that in the case of 
RAYMODE vs PE (Fig. 14) at an FOM of 70 the average difference 
not including the A1SSH run is only 4.0%. At FOM's of 65 and 
70 dB the average difference in detection probability when 
compared to PE is fairly close. At an FOM of 75 aB however, 
ASTRAL vs PE has an average difference of 45.1% while RAYMODE 
vs PE is only 7.9%. The ranges at an FOM of 80 dB were in 
excess of 15 nm and resulted in CDP's of 100 percent in all 
three models with the exception of the cross layer high 
frequency scenarios in the RAYMODE model. (The 80 FOM results 
were not included in CDP comparisons for Area 1.) 

No appreciable difference in computation time between 
ASTRAL and RAYMODE was found using a source frequency of 
either 50 or 400 Hz. Typically both models required 
approximately 3 to 5 seconds to generate proploss curves at 50 
Hz. At 400 Hz the ASTRAL model took an average of 2 seconds 


longer to generate the propagation loss curve than the RAYMODE 


model which generally took 5 sec. 
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Figure 9 - Proploss for Area 1: 400ft/450ft/400 Hz (a and b) 
400ft/450ft/50 Hz (c and d) 
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(dashed) 400 ft/150 ft/400 Hz 


Range (nm) 
RAYMODE (solid) & PE (dashed) 400 ft/150 ft/400 Hz 


Range (nm) 
RAYMODE (solid) & PE (dashed) 400 ft/150 ft/50 Hz 


Range (nm) 


Figure 10 - Proploss for Area 1: 400ft/150ft/400 Hz 
(a and b), 400ft/150ft/50 Hz (c and qd) 
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ASTRAL (solid) & PE (dashed) 90 ft/450 ft/400 Hz 


Range (nm) 
RAYMODE (solid) & PE (dashed) 90 ft/450 ft/400 Hz 


i 
H 
i 


Range (nm) 


ASTRAL (solid) & PE (dashed) 90 ft/450 ft/50 Hz 


Range (nm) 
RAYMODE (solid) & PE (dashed 





Range (nm) 





Figure 11 - Proploss for Area 1: 90ft/450ft/400 Hz 
(a and b), 90ft/450ft/50 Hz (c and d) 
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ASTRAL (solid) & PE (dashed) 90 ft/150 ft/400 Hz 


Range (nm) 
RAYMODE (solid) & PE (dashed) 90 ft/150 ft/400 Hz 


Range (nm) 


Range (nm) 
RAYMODE (solid) & PE (dashed) 90 ft/150 ft/50 Hz 


{ 





Range (nm) 





Figure 12 - Proploss for Area 1: 90ft/150ft/400 Hz 
(a and b), 90ft/150ft/50 Hz (c and d) 
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|__Run__JASTRAL| PE | DIFF [| Run | RAY [PE | DIFF _| 
PAIMMH | 1.90 | 3.412 
PAIMML [4.80 3.447 PAIMML | 2.40 [3.4 [O07 | 
| AIMSH | 48 | 35 [1.3 | |AIMSH| 2.5 | 3.5 [-1.0_ | 
| AiMSE_| 5.1 | 3.4 [2.0 PAIMSL | 2.5 | 3.4 [0.6 | 
|_AISMH | 48 | 25 | 23 | LAISMH | 14.4 | 25 [1.4 | 
PAISML | 1.20 [ 240 | -t2 


A1SSH A1SSH | 2.5 
AiSSL 


| 9.6 | -7.7 
ASSL 


J__Run__JASTRAL| PE | DIFF _| {Run | RAY [PE | DIFF _| 
|_AIMMH | 16.8 | 8.7 [8.1 | [AIMMH | 6.10 | 87 | 3.6 
AIMML | 25.2 | 7.7 {47.5 | 
|_AiMSH | 140 [ 96 | 44 | 
|_AISMH [19.2 | 9.6 | 9.6 | 
| AISML_ | 20.4 | 7.7 [12.7 | 
| AISSH | 75.4 {| 823 | 69 | 
Avg diff | 10.8 _| 


5 

| Run ASTRAL] PE {DIFF _| | Run {| RAY [| PE | DIFF | 
| AIMMH | 65.3 | 17.2 | 48.14 | 
| AIMML | 94.1 | 30.1 | 64.0 | LAIMML | 32.1 | 30.1 [20 | 
| AIMSH | 27.8 | 16.8 | 11.0 | AiMSH 
| 42.1 | 
| AISML_ [| 98.1 | 28.6 |” 69.5 | 
|_AiSSH_ | 100.0 | 100.0 [7 0.0 | 
|_AiSSL_| 100.0 | 20.4 | 79.6 | 

Avg diff | 45.1_| 





Figure 13 - Cumulative Detection Probability Comparisons (%) 
ASTRAL-PE (left) and RAYMODE-PE (right) 








B. AREA 2: GULF OF OMAN 

The transmission losses in area 2 were much greater than 
in Area 1. The propagation loss profiles in Figs. 14-17 show 
a greater than 80dB loss in less than five miles in all but 
the 90ft/150ft/400 Hz case. The rapid loss of sound energy in 
such a short range resulted in a much higher correlation 
between all three models. (The cyclical appearance of the 50 
Hz PE proploss profiles, seen here in Area 2 and even more 
apparent in Area 3, is examined further in Appendix B). The 
average differences in the predicted median detection ranges 
are shown in Table III. ASTRAL consistently predicts longer 
ranges than PE, while RAYMODE generally predicts ranges 
shorter than PE. 

The small differences in ranges similarly resulted in CDP 
differences (Fig. 18) usually less than 10%. Just as in Area 
1 computations, no significant difference in the computation 


time between ASTRAL and RAYMODE was noted. 
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ASTRAL- PE RAYMODE - PE ASTRAL - RAYMODE 


AVG DIF 


| Run | FOM75__| 
PA2MMH J 4.2 
LA2MML | 0.9 
}A2MSH | 4.2 
PA2MSL_ {| 1.0 
LA2SMH | 1.0 
LA2SML | 0.9 
LA2SSL_| 0.9 


AVG DIF 


| Run | _FOM80 
A2MMH | 16. +! 





Table III - Area 2 Average Predicted Detection Range 
Differences 





ASTRAL (solid) & PE (dashed) 400 £t/450 ft/400 Hz (a) 
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Range (nm) 
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ASTRAL (solid) & PE (dashed) 400 ft/450 ft/50 Hz (c) . 
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Figure 14 - Proploss for Area 2: 400ft/450ft/400 Hz 
(a and b), 400ft/450ft/50 Hz (c and 4d) 
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Range (nm) 
RAYMODE (solid) & PE (dashed) 400 ft/150 ft/400 Hz 


Range (nm) 
RAYMODE (solid) & PE (dashed) 400 ft/150 ft/50 Hz 





Figure 15 - Proploss for Area 2: 400ft/150ft/400 Hz 
(a and b), 400ft/150ft/50 Hz (c and qa) 
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ASTRAL (solid) & PE dashed) 90 ft/450 ft/400 Hz 


10 
Range (nm) 
RAYMODE (solid) & PE (dashed) 90 ft/450 ft/400 Hz 
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Figure 16 - Proploss for Area 2: 90ft/450ft/400 Hz 
(a and b), 90ft/450ft/50 Hz (c and qd) 
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ASTRAL (solid) & PE (dashed) 90 ft/150 ft/400 Hz 
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Figure 17 - Proploss for Area 2: 90ft/150ft/400 Hz 
(a and b), 90ft/150ft/50 Hz (c and qd) 
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FOM 65 
| Run ASTRAL] PE | DIFF _| 
FA2MMH [4.32.5 1.8 
A2MML_| 3.4 | 3.5 | -0.4 
A2MSH | 5.4 | 25 | 26 | 
AQMSL_| 4.3 [3.4] 1.2 


| 3.4 | 
| 3.1 | 
| 6.2 | IA2SSH_| 3.5 | 
3.5 | A2SsL_|_7.7 

Avg aift | 


| 4.2 | 
| 0.2 | 
| 1.3 | 
Ltt | 


|A2MSH | 17.0 | 6.2 {| 10.8 | 
PA2MSL | 17.2 | 87 {[ 85 | 
PA2SMH | 17.2 | 8.7 | 8.5 | 
PA2SML {| 17.0 [| 87 | 83 | 
|} A2SSH | _ 67.9 | 74.3 [| -6.4 | 
PA2SSL_|_ 17.2 | 9.6 | 76 _| 





| Run_ ASTRAL} PE | DIFF _| | Run | RAY | PE | 

JA2MMH | 24.8 | 8.7 | 16.1 | [A2MMH| 6.8 | 8.7 

PA2MML | 36.5 | 19.2 | 17.3 _| 

|A2MSH | 42.1 | 16.8 | 25.3 | JA2MSH | 6.8 | 16.8 | -10.0_| 

}A2MSL | 30.4 J 17.2 4} 13.2 | }A2MSL | 6.2 | 17.2 | -11.0 | 

|A2SMH | 49.3 | 17.2 {32.1 _ | |A2SMH | 5.8 | 17.2 | -11.4 | 

| A2SML_{_30.4 | 19.2 | 44.2 | |A2SML | 5.8 {19.2 | -13.4 | 

| A2SSH | 98.6 | 100.0 | -14 | | A2SSH | 100.0 _[ 100.0 | 0.0 | 

LA2SSL | 30.4 17.2 [13.2 _| | A2SSt | 50.0 | 17.2 [32.8 | 
Avg diff {| 16.2 _| Avg diff | 11.4 | 





Figure 18 - Cumulative Detection Probability Comparisons (%) 
ASTRAL-PE (left) and RAYMODE-PE (right) 








C. AREA 3: SOUTH CHINA SEA 

The propagation loss profiles in Area 3 (Figs. 19-22) 
showed very large waeiavionse In the first two scenarios 
ASTRAL and RAYMODE predict ranges two times greater than those 
generated from PE for a FOM of 70 @B. At an FOM of 75 dB the 
ranges are as much as five times greater. Even though both 
models differ greatly from PE the correlation between ASTRAL 
and RAYMODE is better than in areas 1 and 2. -In the 400 Hz 
runs the RAYMODE and ASTRAL profiles have a much higher 
agreement with PE than in the 50 Hz runs. The average 
differences for just the 50 Hz runs for ASTRAL vs PE and 
RAYMODE vs PE, at an FOM of 70, are 2.1 nm and 2.3 nm 
respectively. For the 400 Hz runs the differences are only 
1.1 and .55 nm respectively. Similar relationships are true 
tor the other FOM’s as well. The reason for this is not 
exactly known. The overall differences in ranges are shown in 
Table IV. 

The differences between the detection probabilities (Fig. 
23) for all three models appear to be similar in this area. 
Approximately 2% at FOM of 65 aB, and 13% at FOM of 70 GB. 
The big difference in this area is that correlation between 
ASTRAL and RAYMODE is better overall than either ASTRAL to PE 


or RAYMODE to PE. 
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ASTRAL - PE 


RAYMODE - PE 


FOM 70 


1.7 
1.4 
0.7 


FOM 75 
J 88 
ee 
a 
a: ee 
a) ee 
er eee 
a ae 
a 

FOM 80 


ASTRAL - RAYMODE 


FOM 70 





Table IV - Area 3 Average Predicted Detection Range 


Differences 











Range (nm) 


RAYMODE (solid) & PE (dashed) 400 ft/450 ft/400 Hz 





~ Range (nm) 
RAYMODE (solid) & PE (dashed) 400 ft/450 ft/50 Hz 






Range (nm) 





Figure 19 - Proploss for Area 3: 400ft/450ft/400 Hz 
(a and b), 400ft/450ft/50 Hz (c and da) 
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Range (nm) 


ASTRAL (solid) & PE (dashed) 400 ft/150 ft/50 Hz 


Range (nm) 
RAYMODE (solid) & PE (dashed) 400 ft/150 ft/50 Hz 


Range (nm) 








Figure 290 ~- Proploss for Area 3: 400ft/150ft/400 Hz 
(a and b), 400ft/150ft/50 Hz (c and ad) | 
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ASTRAL (solid) & PE (dashed) 90 ft/450 ft/400 Hz 


10 
Range (nm) 
RAYMODE (solid) & PE (dashed) 90 ft/450 ft/400 Hz 
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ASTRAL (solid) & PE (dashed) 90 ft/450 ft/50 Hz 
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Range (nm) 
RAYMODE (solid) & PE (dashed) 90 ft/450 ft/50 Hz 


Range (nm) 


Figure 21 - Proploss for Area 3: 90ft/450ft/400 Hz 
(a and b), 90ft/450ft/50 Hz (c and ad) 
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ASTRAL (solid) & PE (dashed) 90 ft/150 ft/400 Hz 
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Figure 22 ~ Proploss for Area 3: 90ft/150ft/400 Hz 
(a and b), 90ft/150ft/50 Hz (c and d) 
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[| Run_JASTRAL] PE | DIFF 
|A3MMH | 12.9 | 43 | 8.6 | 
PASMML | 17.2 | 4.8 | 12.4 | 
LA3MSH |_14.0 | 3.5 | 10.5 | 
| ABMSL | 17.2 1 4.30 [12.9 | 
|A3SMH | 16.8 | 3.1 | 13.7_| 
| ASSML_| 17.3 | 3.4 | 14.2 | 
| ASSSH | 23.3 | 248 | -15 | 
LABSSL | 33.8 | 3.1 | 30.7__| 


| Run _JASTRAL] PE | DIFF | 
}A3MMH | 23.3 [5.8 | 17.5 
LASMML | 98.1 | 9.6 | 88.5 _| 
| ASMSH | 33.8 | 62 {| 27.6 | 
| A3SMH |_27.8 | 6.2) | 21.6 | 
}ASSML | 725 | 68 | 65.7 | 
LASSSH |_ 87.8 | 716 | 16.2 | 
LA3SSL_| 93.2 | 6.2 | 87.0 _| 

[Avg diff | 50.2 _| 


| Run _JASTRAL[ PE | DIFF _| 
LA3MML {100.0 [75.4 [| 246 | 
LASMSH | 89.9 {| 55.7 [34.2 | 

| ASMSL | 100.0 | 60.1 _| 
| A3SMH | 87.8 | 73.2_ | 14.6 _| 
| ASSML | 100.0 [| 82.3 [17.7 | 
A3SSH | 97.0 | 984 | -1.1_ | 
A3SSt_| 100.0 | 65.3 | 34.7_| 
Avg diff | 34.9 | 


A3SSL_j| 16.8 | 3.1 | 13.7 _ | 


| Run | RAY | PE | OONFF 
}ASMMH | 57.1 | 5.8 | 51.3 | 
}A3MML {99.7 | 9.6 | 90.1 | 
|ASMSH | 8.2, | 6.2 | 2.0 | 
}ASMSL | 71.6 | 10.7, | 60.9 | 
71. 


| Run | RAY {| PE | DIFF | 
|ASMMH | 94.1 {| 10.9 | 83.2 | 
LA3MML {100.0 | 75.4 | 246 | 
|} A3MSH | 32.1 | 55.7 | -23.6 | 
| A3MSL_| 100.0 | 60.1 [39.9 | 
LA3SML | 100.0 {| 82.3. | 17.7 | 
| A3SSH | 97.8 | 981 | -03 | 
LASSSL | 95.1 | 65.3 | 29.8 | 
Avg diff | 





Figure 23 - Cumulative Detection Probability Comparisons (%) 


ASTRAL-PE (left) and RAYMODE-PE (right) 











V. CONCLUSIONS 

The ASTRAL propagation loss model consistently predicted 
much longer detection ranges than either the PE or RAYMODE 
models. This could lead to incorrectly spaced patterns and an 
overly optimistic probability of detection. RAYMODE 
generally predicted ranges consistent with or slightly shorter 
than PE. This would result in a more conservative search and 
provide a lower chance of allowing a target to slip through an 
area undetected. At lower figures of merit the difference 
between models is fairly small. However, at higher FOM’s 
there are significant differences in transmission loss 
predictions. 

No attempt was made to compare the ASTRAL or RAYMODE 
predictions with empirically gathered data, thus no 
conclusions as to which model better reflects reality is made. 

With the current computer technology (ASQ 212 computer) 
being installed on ASW aircraft, processing capability and 
storage memory for new decision support systems are limited. 
Considering this limitation combined with the evolution of 
quieter submarines and shorter detection ranges, installing 
the more complex and computation intensive ASTRAL program may 
not be the best choice. The passive range independent RAYMODE 
propagation loss model, as shown in this report, will provide 


adequate transmission loss predictions. 
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APPENDIX A 


A. RUN IDENTIFICATION 

Each propagation loss run was identified by a consistent 
code as follows: {A num B C D E}, where A represents Area, 
"num" is either 1, 2, or 3, B is the receiver depth (M=400ft, 
S=90ft), C is the source depth (M=450ft, S=150ft), D is the 
source frequency (H=400 Hz, L=50 Hz), and E is either RD, RI, 
or PE (RD=ASTRAL, RI=RAYMODE, PE=Parabolic Equation Model). 
For example A1MMHRD is a run in Area 1 with receiver at 400ft, 
source at 450ft, frequency of 400 Hz, and was run using 


ASTRAL. 


B. ASWTDA DATA COLLECTION 

Data from the ASWTDA PE and ASTRAL models was accessed 
after the propagation loss plot was generated. The system 
stored the current graph data in a file called OUTPUT.LOG in 
the exec/data/tactical directory for ASTRAL and in a file 
called for042.dat in the exec directory for PE. By opening a 
new window and changing to the necessary directory the files 
were copied to a new file. If this was not done the next 
graph created would overwrite the previous information. The 
ASTRAL data was in ASCII format and immediately available. 
The PE information was in a binary form and had to be 


converted via a Fortran program listed in Appendix B. The 
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RAYMODE data was not stored for retrieval as in ASTRAL or PE. 
Instead, using the "TDA", "Range Independent", "Proploss", 
"Info" selections form the menu bar, a tabular summary of 
RAYMODE range and propagation loss is displayed. By using the 
mouse and left button to highlight the vectors, then "vi" to 
open an editor in an open window, the center mouse button was 
used to copy the data to a new text file with the appropriate 
name. (Ref. 7] | 

The proploss files were edited and made into "m" files 
for the purpose of creating and analyzing the graphs using the 
Matlab program. 

The following tables contain the predicted Median 


detection ranges from all three models at each of the four 


Figures of Merit (65, 70, 75, 80 QB). 








PE 
RUN/FOM [65 [70 | 75 | 80 
PEMMHRI 


| 20.6 | 30.0 | 
PESSLRI 





Table V - Area 1 Median Detection Ranges 
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A2MSLRI_ | 0.7 | 1.0 | 16 | 
A2SMHRI_| 0.7 | 1.0 | 14 | 





Table VI - Area 2 Median Detection Ranges 








Table VII - Area 3 Median Detection Ranges 





APPENDIX 2 


A. PE DATA CONVERSION 

The following Fortran program was used to convert ASWTDA 
binary format (PE output) to ASCII output. Load into any 
directory and run by calling filename.[Ref. 8] 


character*32 infile, outfile 
character*9 when 
integer*2 ititl(40),nr, nz pence: nbotm, nthb,nd 
real £,zs,x%Z(20),xr, thmin, thmax, the,d1,d2 
real tthb(20) ,rthb(20) ,rplot, t1(20) 
print *,’Enter input file (for042.dat format) name:’ 
read(5,*)infile 
open (unit=1, form=’unformatted’ ,status=’old’,file=infile) 
print *,’Enter output file name:’ 
read(5,*)outfile 
open(unit=2,status=’new’ ,file=outfile) 

c READ THE HEADER RECORD OF THE TL FILE 
read(1)ititl,when,f,zs,nr,nz, (x2(i),i=1,nz) 

c write out four lines of header information 
print *,ititl 
write(2,¥*)ititl 
print *,when,f,zs,nr,nz, (xz(i),i=1,nz) 
write(2,*)when,f,zs,nr,nz, (xz(i),i=1,nz) 
print *,’Range TL’ 
write(2,*)’Range TL’ 


print 4, lee meee enn wee ese ren : 
write (2,%) ! q----------- 8 - ee - , 
nd=1 
print *,’Number of depths is’,nd 
c read the range and tl values, print them to screen and c 


to the output file 
do 33 j=1,400 
read(1,end=888,err=999)rplot, (t1(i),i=1,nd) 
write(2,*)rplot,tl(i),i=1,nd) 
33 continue 
888 print *,’End of file on .042 file’ 
close (1) 
close (2) 
print *,’Normal termination’ 
stop 
999 close(1) 











close (2) 

print *, ‘Abnormal termination’ 
stop 

end 


B. NESST 

The program used to evaluated the Cumulative Detection 
Probabilities of the three patterns of interest was the NORDA. 
Environmentally Sensitive Sonobuoy Tactics (NESST) program, 
provided by Wagner and Assoc. NESST has two main 
capabilities. First, the performance of an arbitrary sonobuoy 
pattern can be evaluated with respect to a user specified 
environment. Second, for a given environmental description of 
a region, NESST will develop a sonobuoy pattern which 
optimally searches that region. [Ref. 11] 

Only the first of these two capabilities was used for 
this project. For each of the search regions a 2x8,4x4, and 
5x6x5 sonobuoy pattern was constructed for each FOM, using 1.5 
times the MDR for sonobuoy spacing and 1.67 times MDR for row 
spacing. To evaluate the sonobuoy patterns NESST records the 
location of each buoy in the pattern in terms of grid 
coordinates. The NESST program used the environment created 
by ASWTDA for the 50 nm by 50 nm regions as inesbanis for the 
grid. The patterns were centered on the area of interest and 
evaluated for probability of detection based on the full 50 x 
50 nm area. The sonobuoys were evaluated as cookie cutter 


detectors with the detection ranges predicted by the given 
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model for the grid location. When the full region of interest 
was 100 percent covered by the cookie cutter detectors, the 
pattern was assumed to have a 100 percent probability of 


detection. 


C. CYCLIC EFFECT IN PE PROFILES 

The cyclic appearance seen in the 50 Hz PE proploss 
profiles is due to bottom reflections. The PE model generates 
transmission loss information using a +/- 20° beamwidth. In 
areas where the sound speed profile is strongly negative, like 
in Area 2 and Area 3, this narrow beamwidth is refracted 
downward and reflected off the bottom as seen in Fig. 24. 
When the beamwidth is increased to +/- 40° (see Fig. 25) the 
cyclic effect disappears. Both the ASTRAL and RAYMODE models 
use the wider beamwidth for their calculations. This seems to 
indicate that in cases where there is a large amount of bottom 
interaction the PE model is less than ideal. Since the 
ranges examined in this report are generally less than 10 nn, 


the information presented was not affected by this cyclic 


effect. 








Figure 24 - Range Dependent Raytrace 
with Beamwidth of +/- 20° 
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Figure 25 - Range Dependent Raytrace with 
Beamwidth of +/- 40° 
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